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Abstract. We review recent developments on quantum scattering from 
mcsoscopic systems. Various spatial geometries whose closed analogs shows 
diffusive, localized or critical behavior are considered. These are features that 
cannot be described by the universal Random Matrix Theory results. Instead one 
has to go beyond this approximation and incorporate them in a non-perturbative 
way. Here, we pay particular emphasis to the traces of these non-universal 
characteristics, in the distribution of the Wigner delay times and resonance 
widths. The former quantity captures time dependent aspects of quantum 
scattering while the latter is associated with the poles of the scattering matrix, 
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1. Introduction 

Quantum mechanical scattering in systems with complex internal dynamics has been a 
subject of intensive research activity for a number of years. This interest was motivated 
by various areas of physics, ranging from nuclear £Q, atomic [2] and molecular [3] 
physics, to mesoscopics @], quantum chaos p3E], and classical wave scattering 
Recently, the interest in this subject was renewed due to technological developments 
in quantum optics associated with the construction of new type of lasers Elj and 
experimental investigation of atoms in optical lattices 

The most fundamental object characterizing the process of quantum scattering is 
the unitary S -matrix relating the amplitudes of incoming waves to the amplitudes of 
outgoing waves. The recent advances in mesoscopic physics and quantum chaos, led 
to a fast development of powerful theoretical techniques which allow us to understand 
the statistical properties of the S— matrix. At present, there are two complementary 
theoretical tools employed to calculate statistical properties of the S— matrix, namely 
the semiclassical and the stochastic approach. The starting point of the first is a 
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representation of the S— matrix elements in terms of a sum of classical orbits [5] E] 
while the latter exploits the similarity with ensembles of Random Matrices At the 
same time, recent experimental progress allowed a direct comparison of the theoretical 
predictions with actual experimental results. Microwave experiments (see 114) 
and references therein) offer a unique possibility to check even details of the existing 
theories in cases, where this is hardly possible by other methods while they raise new 
challenging questions (see for example [TTfl 1X41115] 

Apart from the study of the distribution of the S— matrix elements, resonance 
widths and Wigner delay times distributions gained also a lot of attention. The 
latter quantity captures the time- dependent aspects of quantum scattering. It can 
be interpreted as the typical time an almost monochromatic wave packet remains in 
the interaction region. Resonances are defined as poles of the S'-matrix occurring 
at complex energies S n = E n — jT n , where E n is the position and r n the width of 
the resonance. They correspond to " eigenstates" of the open system that decay in 
time due to the coupling to the "outside world" and they are related to conductance 
fluctuations and current relaxation ^B). For chaotic/ballistic systems Random Matrix 
Theory (RMT) is applicable, and the distributions of resonance widths V(T) and 
Wigner delay times V(t) are known. A review of the RMT results can be found in 
PI (see also 0). 

In this contributions we aim in giving an overview of the recent developments 
in scattering from open samples in conditions where RMT is not applicable and 
deviations from "universality" due to the appearance of localization is apparent. 
Although our presentation is focus on random media, one has to keep always in mind 
that these results valid also for dynamical systems with chaotic classical limit. Despite 
the fact that these systems are deterministic (in contrast to random media where 
randomness is "builded" up with the system) localization occurs due to complicated 
interference effects created by the underlying classical chaotic dynamics and for this 
reason is termed dynamical localization |17| . 

The observables that will be in the focus of our presentation are the distributions 
of resonance widths and delay times. We consider various spatial geometries and 
models whose closed analogs show features such as diffusion, criticality or localization. 
A short overview of localization theory and the definitions of the various regimes is 
given in Section [3 In the next Section El we present the mathematical formalism 
associated with the scattering process, and define the quantities of interest. In Section 
0]we review the consequences of localization in the resonance width distribution. The 
corresponding results for the delay times are analyzed in [5J Finally in Section we 
present some results for quasi-periodic systems at criticality. Our conclusions are given 
at the last Section 

2. Various Regimes in Localization Theory: A brief overview 

Localization of waves has always been among the most difficult yet most fascinating 
topics in the study of wave propagation in disordered media. The first studies dealt 
with infinite media, showing that localization is always achieved in one and two 
dimensions but that a minimum amount of disorder is required in dimensions larger 
than two ^3 EH HOI' Its main feature is that the eigenfunctions of a disordered 
medium in the localized regime, are characterized by an exponential decay in space 
i.e. j^r^r)! ~ exp(— |r— ro|/£), where £ is the localization length. A direct consequence 
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is that transmission is inhibited and the system behaves as an insulator. 

The fingerprints of localization in various quantities associated with the closed 
systems have been identified and quite well understood. Detailed numerical and 
theoretical studies gave a clear picture how the statistical properties of these quantities 
change as the disorder strength increase (for a recent review see [21] and references 
therein). For weak disorder, such that the mean free path / mean is larger than any 
physical length scale, the system is in the ballistic regime. An equivalent definition 
of this regime is given through the condition g = £rh/A — > oo. The dimensionless 
ratio g is known as the Thouless conductance, A is the mean level spacing and i/ph 
is the co-called Thouless energy. The latter plays a prominent role in the theory 
of random media, and is inverse proportional to the time needed for an excitation 
to propagate through the entire sample. In this regime the predictions of Random 
Matrix Theory (RMT) were shown to describe accurately the statistical properties of 
various observables. For example, the eigenstates are extended all over the system, 
the eigenvalue spacing distribution follow with a good accuracy the famous Wigner 
surmise [221 [22] etc. 

As the disorder increases the system becomes diffusive. This regime is 
characterized by the condition that the system size L is larger than the mean free 
path / m ean but still smaller than the localization length £ i.e. Z mean « L « ^. Using 
the powerful sigma-model approach it was explained how the deviations from the 
RMT results raise [^ [^] . Detailed numerical experiments [251 1261 [23 verified the 
theoretical predictions. These deviations are specifically strong at the far "tails" of the 
distribution of the eigenfunction intensities as well as of some related quantities and 
are signatures of the underlying classical diffusive dynamics [23 123 12Z] ■ They were 
shown to be related with anomalously localized states termed pre-localized states. In 
[231211 it was found that pre-localized states are also present in quantum systems with 
deterministic chaotic dynamics. As far as the spectral correlations are concerned it 
has been shown that there are large deviations above the Thouless energy Erh = H/td 
where TTh = hD/L 2 is time to diffuse through the system with diffusion constant D 
[212 EEl- The Thouless conductance g is related with the latter as g = DL d ~ 2 where 
d is the dimensionality of the system. Rapid development in microwave techniques 
allowed for a direct observation of some of these predictions, in microwave experiments 

[23113 EH E21 

The deviations of the level and eigenfunction statistics from their RMT form, 
strengthen with increasing disorder and become especially pronounced at the 
localization regime . In this regime, inhibition of wave diffusion due to interference 
of multiple scattering waves take place ^SJ. The resulting scenario depends strongly 
on the dimensionality of the sample and the disorder strength. It turns out that the 
eigenstates are exponentially localized in low dimensional systems even for arbitrary 
weak disorder. As a matter of fact the localization regime is defined through the 
condition that L < £. One of the consequences following from this fact is the prediction 
that the conductance of a sample goes exponentially to zero with increasing its length 
and the sample behaves as an insulator. In contrast, disordered single-particle systems 
in more than two-dimensions exhibits a readier behavior. If the disorder is weak 
enough there is no localization and the system has a metallic behavior while for 
strong disorder strength we recover the localization regime. The transition point from 
a metallic to localized behavior is of special interest and is called Metal-Insulator 
Transition (MIT). 

The MIT where the phase transition from localized to extended states occurs, is 
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characterized by remarkably rich critical properties. In particular, the level statistics 
acquires a scale-independent form with distinct critical features [21 HHI IHH] while the 
eigenfunctions show strong fluctuations on all length scales and represent multi-fractal 
distributions [211 1241 l3"H 1351 136) . As a matter of fact, a connection between multi- 
fractality and statistical properties of eigenvalues at MIT has been recently established 
|33| . The multi- fractal structure of the eigenfunctions is usually quantified by studying 
the size dependence of the so-called participation numbers (PN) 

Af q = (^J |"0(r)| 29 drj aclM-^ 13 * (1) 

where L is the linear size of the system and D q are the multi-fractal dimensions of the 
eigenfunction ip(r). Among all the dimensions, the correlation dimension D2 plays the 
most prominent role. The corresponding PN is roughly equal to the number of non- 
zero eigenfunction components, and therefore is a good and widely accepted measure 
of the extension of the states. At the same time, D2 manifest itself in a variety of 
other physical observables. As examples we mention the statistical properties of the 
spectrum |21l I33| . and the anomalous spreading of a wave-packet, and the spatial 
dispersion of the diffusion coefficient [2Zj ■ 

At the same time a considerable effort was made to understand the shape of the 
conductance distribution V{g) at MIT 16, 38, 39 . However, it is still unclear whether 
the limiting V (g) is entirely universal, i.e. dependent only on the dimensionality and 
symmetry class, as required by the one-parameter scaling theory of localization [T%] . 
The latter is one of the major achievements in the long history of studying the MIT. 
Its basic assumption is that close to the MIT the change of the conductance g with the 
sample size L depends only on the conductance itself, and not separately on energy, 
disorder, size and shape of the sample, the mean free path etc. 

Although a lot of studies have been devoted to the analysis of eigenfunctions and 
eigenvalues and of conductance, the properties of resonances, and delay times were 
left unexplored until recently. Nevertheless it was clear from the very beginning that 
their statistical properties depend strongly on the nature of the states of the finite 
system " in isolation" . Thus Anderson localization must leave its fingerprints in these 
quantities which after all reflect the "leakage" of the waves to the leads, through the 
sample boundaries. In the next sections we will review the outcome of these studies 
and their deviations from the RMT predictions due to non-universal features. 



3. Quantum Scattering: Basic Concepts 

The scattering S— matrix relates the outgoing wave amplitudes to the incoming wave 
amplitudes. Assuming M open channels, one can show that the M x M scattering 
matrix can be written in the form 

S(£) = 1 - 2^7^y t - — V; H eS = H - inVV^ (2) 

t - fleff 

Here, H stands for an A— dimensional self-adjoint Hamiltonian describing the closed 
counterpart of the system under consideration, £ stands for the energy of the incoming 
waves, and V is an M x N operator that contains matrix elements coupling the internal 
motion to one of the open M channels. In principle, the matrix elements of the 
operator V depends on energy. However, since this dependence is very weak (far away 
from channel thresholds), we can ignore it and consider Vij to be energy-independent. 



Statistics of Resonances and Delay Times in Random Media: Beyond Random Matrix Theoryb 



1 is the MxM unit matrix. For a detail derivation of the scattering matrix for the case 
of a tight-binding model see 0U] while for a chaotic cavity see 0T] . It is easy to verify 
that form @ ensures the unitarity of the scattering matrix i.e. S — 1, provided 
the energy £ takes only real values. When one allows the energy parameter to have a 
nonzero imaginary part, the S— matrix unitarity is immediately lost. Having Xm£ > 
correspond to the physical situation of uniform dumping inside the system E2] and 
it is responsible for the losses of the outgoing flux of the particles as compared to 
the incoming flux. The "dual" case Xm£ < correspond to uniform amplification. 
The balance between the two fluxes is precisely the physical mechanism behind the 
S— matrix unitarity. 

The poles of the scattering matrix S are associated with the formation of 
resonance states. They represent long-lived intermediate states to which bound states 
of a closed system are converted due to coupling to continua. Due to causality, they are 
located in the lower half plane i.e. £ n — E n — \T n , where E n and T n are the position 
and the width of the resonances, respectively. They are solutions of the following 
secular equation 

det(f - H cS (£ )) = 0. (3) 

where the resonance width T is inverse proportional to the lifetime of the corresponding 
resonance state. From Eqs. (12IMI) it is clear that the formation of resonances is closely 
related to the internal dynamics inside the scattering region which is governed by Hq. 

Another quantity that will be in the focus of this contribution is the Wigner delay 
time |43| and its variations (for an overview on the various definitions of delay times 
and their physical importance see 01]). It captures the time-dependent aspects of 
quantum scattering. It can be interpreted as the typical time an almost monochromatic 
wave packet remains in the interaction region. Formally the Wigner delay time tw is 
related with the energetic derivative of the total phase of the unimodular S— matrix 

tw(E) = ^TrQ(£0; Q(E) = -ihS^E)^- (4) 

where Q(E) is called the Wigner-Smith time delay matrix 03]- Its eigenvalues r q , 
are called proper delay times, and correspond to the time the particle dwells at a 
particular channel q = 1, • • • , M. Alternatively, one can also define the partial delay 
times tP as the energetic derivatives of the eigenphases {9 q },q = 1,---,M of the 
unimodular S— matrix i.e. tP = d6 q /dE Beyond the one-channel case, proper 

and partial delay times differ, although the sum of partial/proper delay times over all 
M scattering channels are always equal and yield the Wigner delay time. 

4. Resonances 

In this section we analyze the distribution of resonances V(T). The properties of 
resonances are of fundamental as well as technological interest. One can show that they 
determine the conductance fluctuations of a quantum dot in the Coulomb blockade 
regime |45| . or the current relaxation. The latter study constitute a fundamental 
source of physical information for systems which are coupled to a continuum via 
metallic leads or absorbing boundaries. While the radioactive decay is a prominent 
paradigm, more recent examples include atoms in optically generated lattices and 
billiards the ionization of molecular Rydberg states 02], photoluminescence 
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spectroscopy of excitation relaxation in semiconductor quantum dots and wires |49|. 
and pulse propagation studies with electromagnetic waves |5(J) . 

From the theoretical point of view, one can approach the problem of current 
relaxation by evaluating the survival probability P(t) of a wave packet which is initially 
localized inside an open sample of volume f2 

p (t)= [ |*(i,r)| 2 dr « ( dTTV (T) exp(-Tt) . (5) 
Jn Jo 

The approximation above (modal approximation) valid I5UI 1511 1521 15*5] for times 
larger that the Heisenberg time and allow us to calculate a dynamical quantity such 
as P(t), based on information about resonances. The total current leaking out of the 
sample is then related to the survival probability by 

The ability of constructing micro-lasers with chaotic resonators which produce 
high-power directional emission [5] as much as the experimental realizations of the so- 
called random lasers ^U] where the feedback is due to multiple scattering within the 
medium (instead of being due to mirrors) is another reason why statistical properties 
of resonances became of fashion in our days. For the latter application the knowledge 
of resonance width distribution can result in the knowledge of the statistical properties 
of the lasing threshold. 

The lasing threshold is given by the value of the smallest decay rate (i.e. 
smallest resonance width) of all eigenmodes in the amplification window [541 155) . The 
underlying reasoning is that in the mode with the smallest decay rate the photons are 
created faster by amplification than they can leave (decay) the sample. Assuming that 
the number of modes K ^> 1 that lies in the frequency window where the amplification 
is possible, have resonance widths V that are statistically independent one gets for 
distribution of lasing thresholds V(T) |54l 1531 : 

-, K-l 



V{T) = KV(T « 1) 



1 - / V(V « l)dV 
Jo 



(7) 



where we have assumed that all K resonances are distributed according to "P(r << 1). 
The validity of this approximation was verified recently in the framework of the RMT 
|57|. An important outcome of |25ll2(:ij was that one can identify in diffusive systems, 
traces of pre- localized states in the latter distribution and consequently in V(T). This 
send some light to recent experimental finding for random lasers which suggests the 
appearance of localized modes in diffusive samples 58 . 



4-1- Ballistic Regime 

At the ballistic regime, RMT modeling is applicable. Its main advantage is its 
universality. At the same time, universality means that the RMT modeling doesn't 
"know" anything about the specific properties of the system under study. Since 
no physical parameters (except the global symmetries) are plugged in the RMT 
machinery, it is clear that it can give the correct predictions only in some limiting 
cases, when all physical parameters and scales can be considered as irrelevant. 
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Figure 1. Distributions of the resonance widths in the ballistic regime for 
two different models: (a) The Open Kicked Rotor |26| and (b) a fully connected 
quantum graph [H] with "generic" vertex scattering matrices. The numerical data 
(o) are in excellent agreement with the RMT predictions of Eq. JSJ (dashed lines). 



In the general case, Fyodorov and Sommers proved that the distribution of 
scaled resonance widths V — T/A for the unitary random matrix ensemble, is given 

by 

where the parameter q controls the degree of coupling with the channels, (• • •) indicates 
an average over realizations and "y(.) is the 7— function. In Fig. ^ we report some 
representative results from two models in the ballistic regime together with the 
theoretical prediction (JSJ. The excellent agreement is evident. 

In the limit ofM> 1, Eq. JSJ) reduces to the following expression ^2] 

M for _M < f < M 



V(T) = < ^ 2 ' 7r (?+ 1 ) . (9) 

[ , otherwise 

The following argument provide some intuition about the form of resonance width 
distribution . First we need to recall that the inverse of T represents the quantum 
lifetime of a particle in the corresponding resonant state escaping into the leads. 
Moreover we assume that the particles are uniformly distributed inside the sample 
and spread ballistically until they reach the boundary, where they are absorbed. Then 
we can associate the corresponding lifetimes with the time t^ ~ l/r^ a particle needs 
to reach the boundaries, when starting a distance R away. This classical picture can 
be justified for all states with T > T c \ ^S> A where r c i is the classical decay rate which 
can be calculated numerically from the exponential decay of the classical probability 
to stay inside the sample. (For RMT models we have the so-called Moldauer-Simonius 
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relation T c i ~ AMln(l — K-S)! 2 ) which give us the lower bound in Eq. © while for 
generic chaotic system T c i ~ H(s/D,)v where s is the width of the opening, f2 the total 
volume of the system and v the velocity of the particle moving inside the system). The 
relative number of states that require a time t < t^ in order to reach the boundaries 
(or equivalently the number of states with T > Tr) of a d— dimensional system with 
linear dimension L is 

*WCT*) - £ v(r )d r ~ ^ = Ld -^- R ) d (10) 

where V(tp) ~ L d — (L — R) d is the volume populated by all particles with lifetimes 
t < tjj. Assuming ballistic motion i.e. R — vtu, we get from Eq. (|10J) in the limit 
where Y > T cl 

Pint(r R ) ~ J- (11) 
1 R 

which eventually leads to the RMT prediction Eq. ©. 

Eqs. (|8I9I) is our starting point. In the next subsections we will investigate how 
deviations from these expressions raise as we increase the randomness of the system. 



4-2. Diffusive Regime 

We start our presentation with the study of small-resonance width distribution 
■p(r < A). The small resonances T < A can be associated, with the existence 
of pre-localized states of the closed system (for a discussion on pre-localized states 
see ESI 1251 121] ) ■ They consist of a short-scale bump (where most of the norm 
is concentrated) and they decay rapidly in a power law fashion from the center of 
localization |21U24ll2l)ll27| . One then expects that states of this type with localization 
centers at the bulk of the sample are affected very weakly by the opening of the system 
at the boundaries. In first order perturbation theory, considering the opening as a 
small perturbation we obtain |261 

~ = = I*HI 2 ~ L d -^(L)\ 2 (12) 

n£boundary 

where ^(L)! 2 is the wavefunction intensity of a pre-localized state at the boundary 
and d is the dimensionality of the sample. At the same time the distribution of 
wavefunction components at the boundary was found to be |24j 

V(9) - cxp (-A ( p ] ln d (0 4 - rf )) , = l/i (d - 1)/2 *(L) (13) 

where the coefficient A^ d ' oc (3D. Here (i = 1(2) denotes the corresponding ensemble 
for preserved (broken) time-reversal symmetry. Using Eq. (|13|l together with Eq. (|12H 
we obtain 

V(T < A) - exp (-Cj d) ln d (l/r)) , where C p cx (3D (14) 

A detailed numerical analysis performed in JJBJ [59J [HUj for the d — 2,3 showed that 
the above perturbative derivation valids as well for relatively large resonances i.e. 
A < T < T c ; where r c / = D/L 2 is the inverse Thouless time. In Fig. [21 we report 
some numerical data for the case d = 2 and compare with the theoretical predictions 
of Eq. (|14|) . Let us finally compare Eq. (|14|) with the results for ballistic systems 
(see Eq. ©). In the latter case, a strip free of resonances is formed. In the diffusive 
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Figure 2. (a) The distribution of resonance widths (plotted as "P{l/T) vs. 1/r) 
for r < r c ; for two representative values of D for the 2- dimensional Kicked 
Rotator in the diffusive regime 

|2I3 GH3- The system size in all cases is L = 80. 
Filled symbols correspond to broken TRS. The solid lines are the best fit of 
Eq. 1141 for /3 = 1(2) to the numerical data, (b) Coefficients Cp vs. D for the 
same model as in (a) EHIESI- The solid lines are the best fits to Cg = ApD + Bp 
for = 1(2). The ratio R = A 2 /Ai = 1.95 ± 0.03. 

regime, in contrast, there are pre-localized states, which are weakly coupled to the 
leads. Due to their existence the distribution of the small resonance widths has a 
non-trivial behavior described by Eq. i|14|) . 

Next we turn to the analysis of V(T) for T > T c /. Using the same argument that 
led to Eq. I|10|) . but assuming now diffusive spreading i.e. R 2 = D x t, we get (to 
leading order approximation with respect to r c i/r) 

V(T) ~ (~) 2 . (15) 

valid for quasi-one [HI], two [251159} or three-dimensional [HOI random media as long as 
the leads are attached to the boundary of the sample. We conclude that the different 
power law decay of Eq. lfl5|l with respect to Eqs. (|9I11|I is a result of the different 
nature of the dynamics: ballistic vs. diffusive. 
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Here it is interesting to point that a different way of opening the system might 
lead to a different power law behavior for V(T). Such a situation can be realized if 
instead of opening the system at the boundaries we introduce "one-site" absorber (or 
one "lead") somewhere inside the sample. In such a case for d = 2 we have 



V(t R ) _ R 2 _ Dtft r. 



7Wr*)~^ = ^ = ^~^ (i6) 



leading to the following power law decay 



P(r)~(i) 2 . (17) 

Similarly, the analog of Eq. 1|17[) in d = 3 is 

P(r)~(i) 2 - 5 (18) 

where we had to substitute V(t R ) ~ R 3 . The above results valid for any number M 
of "leads" such that the ratio M / L d scales as 1/L d . 

If on the other hand we attach the open channel to the boundary (assume square 
geometry) of a 3d sample we come out with a decay law which is the same as the one 
given by Eq. (|17fl . This is due to the fact that the decay from the surface leads to a 
situation alike the one of a 2d system. 

In Fig. [3] we present some numerical calculations of the 2c? Kicked Rotator [5SJ 
while in Fig.^we present numerical data from the 3d Anderson model in the diffusive 
regime |60| . In all cases a comparison with the corresponding theoretical predictions 
I|15I17I18JI shows a nice agreement. 



4-3. Localized Regime 

Various groups |601 1621 1631 IM) had investigated the resonance width distribution in 
the localized regime during the lasts years. 

In the region of exponentially narrow resonances T < T = exp(— 2L/£) the 
distribution was found to be log-normal i.e. 

P(f)^exph(4|)- 1 ln 2 (f)], r<r . (19) 

This is entirely analogous to the conductance distribution of localized systems. 
Equation (|19|) essentially relies on two assumptions: first, that eigenfunction 
components are randomly distributed with no long-range correlations, and second, 
that they are exponentially localized with a normal distribution of localization lengths. 
This part of the distribution becomes negligible at large L, because it comes of a 
fraction ~ £/L of the full set of all resonances. 
Instead the long resonance tails behave as 

^( f )~(|)p r <r«i/L (20) 

Eq. I|20|l can be easily understood once employing Eq. 1)10(1 . The new ingredient now is 
that wavefunctions are exponentially localized i.e. |W(r)| ~ (l/£ d ^ 2 ) exp(— r/£). Using 
simple perturbation arguments, we have that (see Eq. (|12|0 V ~ ^(r)! 2 which leads to 
the following approximation about the volume V(tn) oc R d ~ £ d ln d (£ d r). Inserting 
this in Eq. I|10(l we get (to leading order approximation with respect to £/L) Eq. i|20|) . 




Figure 3. (a) The resonance width distribution V{T) for preserved TRS and 
D = 20.3 (o) and D = 33.5 (o). The corresponding filled symbols represent V(V) 
for broken TRS and the same values of D. The dashed (solid) vertical line mark 
the classical decay rate r c j for D = 20.3(D = 33.5). (b) The "Pi n t(r) for a sample 
with nine leads (lower curve). For comparison we plot also the "Pi n t(r) for the 
same sample but when we open the system from the boundaries. The dashed lines 
correspond to the theoretical predictions 1151 and 1171 . The figure is taken from 

ES3 



The large T region is essentially determined by the coupling to continuum, so it 
should be model-dependent. Nevertheless, it is reasonable to assume that the number 
of resonances involved is constant, of order £ and therefore this extreme tail should 
subside at large L, at a rate ~ £,/L- 

From Eqs. i|19l20fl it becomes evident that T^fT) depends on one parameter. 
Namely the dimensionless parameter This dimensionless parameter is the 

cornerstone of the one-parameter scaling theory of localization. It was shown in the 
past that the dimensionless conductance g is a simple function of £,/L. The above 
theoretical considerations were tested |HS1 EH f° r various disordered and chaotic 
systems with dynamical localizations and were found to describe nicely the numerical 
data. In Fig. we report some representative cases from the 3D Anderson model 60 
in the localized regime. 

Let us finally note that in the thermodynamic limit L — ► oo the probability of 
finding an eigenstate at any finite distance from the boundary is equal to zero. Thus 
the distribution of the resonance widths in this case collapses into a delta function 




Figure 4. (a) The resonance width distribution V{Y) for the 3d Anderson model 
IbUI and various configurations of the open channels. The dashed lines are the 
corresponding theoretical predictions given by Eqs. I15I17I18I . (see text). The 
figures is taken from 1601 . 



centered at zero. 



4-4- Criticality 

The statistical properties of various observables at the MIT is one of the most 
intriguing problems for many years now. Actually, despite the rich activity 
|18l 1211 1241 1281 1341 1361 13"5| very few theoretical results are known. Here, we present 
consequences of the MIT on the statistical properties of the rescaled resonance widths 
r. It was found |65| that V(T) follow a new universal distribution, i.e. independent of 
the microscopic details of the random potential, and number of channels M as can be 
seen from Fig. [5] For small resonance widths i.e. T < 1, we have found (SOj that V(T) 
can be fitted with a log-normal. The sharp peak on the extreme right corresponding 
to very large resonances is non-universal (model specific) and statistically insignificant 
since it subside as L increases like M/L 3 ~ L" 1 (see also discussion at the previous 
section). 

On the other hand, the main part of V(T), corresponding to intermediate large 
resonances, follow a power-law which is different from those found for ballistic, diffusive 
or localized systems (see previous sections) i.e. 

V(f) ^ g 1 J d T- (1+1 / d '>. (21) 

where g c is the Thouless conductance at criticality 

One can relate the power-law decay l|21[) to the anomalous diffusion at the MIT. 
Indeed, at MIT the conductance of a d-dimensional disordered sample has a finite 
value g c ~ 1. Approaching the MIT from the metallic side one has g ~ Ex /A, where 
Et = D/R 2 is the Thouless energy, D is the diffusion coefficient, and A ~ l/R d 
is the mean level spacing in a <i-dimensional sample with linear size R. This yields 
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Figure 5. P(f ) in the localized regime for various combinations of V and L in the 
range T < 1. The log- normal decay is highlighted by Gaussian fits (full curves) 
whose maximum decreases increasing strength of disorder and also shifts towards 
smaller values of T. Keeping the ratio £/L ?s 0.136 fixed, coinciding distributions 
(filled circles and open squares) are obtained for different combinations of L, V. 
(b) For f > 1 the anticipated power-law decay "P(r) ~ 1/f is observed (dashed 
line) which becomes more robust for increasing strength of disorder. The figure 
is taken from 60 . 



D ~ g c /R d 2 at W c . Taking into account that D = R 2 /t, we get for the spreading of 
an excitation at the MIT 

R d (t) ~ g c t . (22) 

A straightforward application of Eq. I|10l) then, leads to Eq. I|21() . In Fig. [7|we report 
some numerical results for the 3D Anderson model at MIT |fi0l I65| . An inverse 
power law Vi n t(T) ~ T~ a is evident. The best fit to the numerical data yields 
a = 0.333 ± 0.005 in accordance with Eq. (Ell- 
in the original proposal of the scaling theory of localization, the conductance g 
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Figure 6. Universal behavior of at the MIT for a 3d Anderson model. The 

figure is taken from 1651 




Figure 7. The integrated distribution Pi n t(f) for the 3d Anderson model at 
MIT 1051101)1 . The dashed line is the theoretical prediction H21I corresponding to 
Pint(f ) ~ f -0 ' 333 for our case. The figure is taken from In0ll65) . 
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Figure 8. (a) T'-^CW, L) as a function of disorder strength W for different system 
sizes L provides a means to determine the critical point W c of the MIT (vertical 
line at W c = 16.5). (b) The one-parameter scaling of V° nt (W,L) [Eq. (23] is 
confirmed for various system sizes L and disorder strengths W using the box 
distribution. The figure is taken from 16011651 . 



is the relevant parameter JHj ■ A manifestation of this statement is seen in Eq. (|21|) 
where = PintC^o) is proportional to the conductance g. It is therefore natural 
to expect that P® nt will follow a scaling behavior for finite L (and for some To ~ 1), 
that is similar to the one obeyed by the conductance g. It was therefore postulated in 
|rjUl IrjS] the following scaling hypothesis 

V& t (W,L)=f(L/Z{W)), (23) 

where £(W) is the correlation length at MIT. In the insulating phase (W > W c ) the 
conductance of a sample with length L behaves as g(L) ~ exp(— £/£) due to the 
exponential localization of the eigenstates, and therefore we have g(L\) < 3(^2) for 
L\ > Lj2- Based on Eq. (|21|1 we expect the same behavior for V® nt i.e. for every 
finite L x > L 2 we must have £>P t (W,Li) < V° nt (W,L 2 ). On the other hand, in the 
metallic regime (W < W c ) we have that g(L) — DL d ~ 2 and therefore for d > 2, we 
expect from Eq. V^ t (W,Li) > T>f nt (W,L 2 ). Thus, the critical point is the one 
at which the size effect changes its sign, or in other words, the point where all curves 
"P^ t (W, L) for various L cross. One can reformulate the last statement by saying that 
in the thermodynamic limit L — ► 00 at W — W c the number of resonances with width 
larger than the mean level spacing goes to a constant. 

In Fig. we show the evolution of Vi nt (W) for different L using the box 
distribution |601 . From this analysis the critical disordered strength W = W c = 
16.5 ± 0.5 was determined in [HOI EE] in agreement with other calculations |18l I3H| . A 
further verification of the scaling hypothesis H23J1 is shown in Fig. |Ht> where the same 
data are reported as a function of the scaling ratio L/£. All points collapse on two 
separate branches for W < W c and W > W c . 
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5. Delay times 

We turn now to the analysis of Wigner and proper delay times as defined in Eq. l|H 
above. Their knowledge is relevant for experiments on frequency and parameter- 
dependent transmission through chaotic microwave cavities jS] I31| or semiconductor 
quantum dots with ballistic point contacts [HHj. Also, it can be shown that they are 
related to the distribution of reflection coefficients R in the present of weak absorption. 

Absorption is one of the main ingredients in actual experimental situations and 
gain a lot of interest the last years (see 02] an d EH in this volume). Unfortunately 
a comprehensive treatment of absorption is still lacking. There are only very few 
reported analytical results for the distribution V(R) of the reflection coefficient 
R = S'S in the presence of absorption and all of them are within the regime of 
applicability of RMT 021 E7J EE], except the recent work EHUZOj, where quasi- ID 
geometry in the localized regime is considered as well. 

Specifically in the weak absorption limit it was shown [HI E7J that the following 
relation holds: 

Rg = 1 - T q /T a (24) 

where r q are the proper delay times (eigenvalues of the Wigner-Smith operator) and 
l/r a is the absorption rate. Thus the knowledge of V(R) reduces to the calculation 
of the distribution of proper delay times V{r q ) 67 . 

For one channel the distribution of delay times is now quite well understood and 
studied in all regimes. Quite recently it was shown |71| the existence of a very general 
relation between the delay time distribution and the distribution of eigenfunction 
intensities 

<^- fc ) = (y k+1 ) (25) 

where y = f2|^> n (r)| 2 is the local eigenfunction intensity and f = tA/2it. Eq. (|25ll 
leads to the following functional relation between the two distributions: 

V w {f w ) = ^V v ( J- J . (26) 

On one hand, this relation allows us to use the existent knowledge on eigenfunction 
statistics |21| to provide explicit expressions for delay times distributions in various 
regimes of interest. On the other hand, since phase shifts and delay times 
are experimentally measurable quantities, especially in the one-channel reflection 
experiment E7J [72] [7JD IZH i this relation opens a new possibility for experimental 
study of eigenfunctions. 

5.1. Ballistic Regime 

We start again our presentation from the ballistic regime. The notion of proper delay 
times goes back more than 40 years to the seminar paper of Smith 031 • Although 
many authors have worked on this problem |12l 1751 1761 1771 1781 17^] , only recently its 
probability distribution was calculated. It was shown in [HO] using standard RMT 
methods, that the distribution of inverse proper delay times is given by the Laguerre 
ensemble from random matrix theory 

V(rr\- ■ ■ , r M *) oc J] K X -rT X f II l(r-y M/2 exp-^^A .( 27 ) 

i<j k 
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As a matter of fact from Eq. i|27|) one can evaluate the distribution of Wigner 
delay times Typ, which for large values decays as a power-law 

1 



V{t w ) oc 



2+/3M/2 



(28) 



in agreement with an earlier conjecture by Fyodorov and Sommers |12|. 

Specifically for M = 1 the distribution of Wigner delay times was calculated even 
in the crossover regime between unitary (/? = 2) and orthogonal (/3 = 1) symmetry 
classes and was found to be [SJJ 



' [fw) "2f3 



d\2 Ao6 



-x\xl-i) 



T 2 (A, A 2 ) = 2X 2 [(1 - X 2 )e- a + Al(l - e 



"(1-0, 



A 2 ^Al^T 



where a = JT 2 (1 — A 2 ), Io(z) stands for the modified Bessel function, and X is a 
crossover driving parameter. For pure symmetries Eq. Ij29(l leads to 1711 ITH] 



P(tw) = [(/3/2)^/ 2 /T(/5/2)]-, 



-/3/2-2 



exp 



-j8/2rw 



(31) 



The following simple argument can be used in order to understand the tails of 
the Wigner delay-time distribution Eq. Our starting point is the well known 

relation 



t(E) 



M 

E 



E n 



r 2 /4 



(32) 



which connects the Wigner delay times and the poles of the 5— matrix. 

It is evident that large times t(E) ~ T^ 1 corresponds to the cases when E ~ E n 
and r n -C 1. In the neighborhood of these points, t{E) can be approximated by a 
single Lorentzian 1|32|) . Sampling the energies with step Ai£ -C T m i n we calculate 
the number of points for which the time delay is larger than some fixed value r (see 
Fig. 13. Assuming that the contribution of each Lorentzian is proportional to its width 
one can estimate this number as <x y T T n /AE. For the integrated distribution of 
delay times in the limit AE ->flwe obtain 



Pint to 



1/7 



dr-p(r)r 



(33) 



and by substituting the small resonance width asymptotic given by Eq. JBJ) we come 
out with the power-law expression l|28l) . 



(30) 



5.2. Diffusive regime 

Using the general relation (|26l) for the case of M = 1, we find for V(t\y) |7T] 

e -i/2f 



V w (t w ) 

7~*w (j~w ) 



1 



4e -2/T 



k / 3 3 1 
4 1 



k / 6 2 , 
2 3 -^ + ^' 



/3 = 1, 
/3 = 2, 
= 4, 



(34) 
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Figure 9. Schematic plot for the Wigner delay time as a function of energy 
according to Eq. 1321 , 



where the parameter k oc g^ 1 is inversely proportional to the dimensionlcss 
conductance g. The proportionality coefficient depends essentially on the sample 
geometry and on the coordinates of the lead. Note that in the limit of g — ► oo we 
recover the RMT results discussed in the previous section. 

Eq. pijl holds for relatively large delay times f w > y^, while in the opposite case 
the distribution is dominated by the existence of the anomalously localized states and 
has the following behavior for dimensionality d — 2,3 |71| : 



1 



< T... < 



(35) 



V w (f w ) ~ exp(-C rf ln d (l/f w )), f w <K. (36) 

Note that although in [7T] the coefhcient Cd was claimed to depend only on the 
dimensionality of the system, it is possible to depend also on the symmetry parameter 
(3 as well (for a discussion on these issues see [21 and references therein). 

For many open channels M ^> 1 there are no quantitative theoretical results 
yet. However one can employ qualitative arguments which together with numerical 
findings can allow us to understand the resulting distributions. Indeed, substituting 
in Eq. I)33|l the small resonance width asymptotic for the V(T) given by Eq. (|14Jl we 
come out with the following log-normal law for the large r regime 

V(T > r d') ~ exp(-a ln d r) (37) 
where the coefficient Cp is the same as the one given in Eq. (|14p. This prediction 
has been tested numerically in 2"> 2(i for the case of a 2d diffusive system, and the 
numerical findings (see Fig. I10f) were shown to be in excellent agreement. 
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Now we estimate the behavior of V(t) for r < r^ 1 . In this regime many short- 
living resonances contribute to the sum l|32[) . We may therefore consider r as a sum 
of many independent positive random variables each of the type t„ = T n x n , where 
x n — 6E~ 2 . Assuming further that 5E n are uniformly distributed random numbers 
we find that the distribution V(x n ) has the asymptotic power law behavior \jx\l" 1 . As 
a next step we find that the distribution V{r n ) decays asymptotically as where 
we use that V(T n ) ~ 1/r^ 2 . Then the corresponding V(t) is known to be a stable 
asymmetric Levy distribution L^\(t) of index // = 1/2 [H2j which has the following 
form at the origin 

Tir^r^^^expi-a/r), (38) 

where a is some constant of order unity. Simple theoretical arguments "T2] suggest 
that this part of the distribution of the Wigner delay times is the same as in RMT 
considerations. This is in contrast with the large delay times (see Eq. (1371) ) where 
RMT considerations lead to a power law decay (J2BJ- 

Since r = Yli=i T qi we expect the behavior of the distribution of proper delay 
times V{T q ) to be similar to V{t) for large values of the arguments (for r>lwe have 
t ~ T™ ax ). The above predictions were verified numerically for a 2d diffusive system 
|2fi| . while unpublished results ["""13] show that the same law applies for a 3d diffusive 
system as well. We point out here that the asymptotic behavior V{t) ~ 1 /r 3 / 2 emerges 
also for chaotic/ballistic systems where the assumption of uniformly distributed 5E n 
is the only crucial ingredient (see for example |12|). 



5.3. Localized Regime 

In a serious of recent works [HH EH ESI EE] it was found that for ID systems with 
M = 1 and weak disorder the delay time distribution is 

V(t w ) = -^exp(-£/vT W ), (39) 

where £ is the localization length and v = \dE/dk\ is the group velocity. (|59")l takes 
its maximum value at r^ ax — Q.5£/v, indicating that the most probable trajectory 
that a particle travels (forth and back) before it scatters outside the sample is the 
mean free path ' mcan = £/4. As tw — > oo, V(tw) shows a long time tail which goes 
as 2t^ clx /t^ v . This leads to a logarithmic divergence of the average value of tw, 
indicating the possibility of the particle traversing the infinite sample before being 
totally reflected. As was indicated in (841 1851 this is another manifestation of the 
fact that in the localized regime the conductance shows log-normal distribution due 
to the presence of Azbel resonances. 

We note that although the distribution for small delay times depends on disorder 
strength and possibly on number of channels M, the long time tail is universal. As 
a matter of fact one can understand the long time power law behavior by employing 
the argument leading to Eq. 133fl . Indeed by substituting the resonance distribution 
in the localized regime 1 (2*0)1 we get again V{t w ) cx 1/t^ independent of the number 
of channels M. The validity of this calculation was checked recently "nUj for the 3d 
Anderson model in the localized regime and for M»l channels. 

The above expression Eq. I)39|l does not contain the length of the chain, indicating 
that this intermediate asymptotics of the delay-time distribution, is related to the 
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Figure 10. The proper delay times distribution "P(r 9 ) for a 2a! Kicked Rotator 
|26l with diffusive coefficients D = 20.3 (o) and D = 29.8 (□ ). The (•) correspond 
to D = 20.3 but now with broken TRS. The dashed lines have slopes equal to C'p 
extracted from the corresponding V(T) (see Fig. 8). In the inset we report Vijq) 
for moderate values of r q in a double logarithmic scale. The figure is taken from 
I2T1I551 . 



resonance width distribution, which is dominated by the electron escape rate from 
the resonant state into the nearest reservoir and for L — > oo is exact for any delay 
time. However, the finite length L determines a cutoff tw ~ exp L ^ for this universal 
behavior, and for larger delay times we find |85| 

T(t w ) oc exp(-L/0r w {1+i lnTw) , r w > exp L /« . (40) 



5.4- Criticality 

Recently, an intensive activity to understand V(tw) for systems at critical conditions 
was undertaken in jSS] 1691 1711 I87| . The activity was mainly concentrated in the 
simplest scattering set up of one open channel attached to the system of linear size 
L. As a result, an anomalous scaling, of inverse moments of Tw with the system 
size L was reported Ej and specific predictions linking the scaling exponents and 
the multi-fractal properties of eigenfunctions of the corresponding closed system were 
made. Specifically it was found that the inverse moments of Wigner delay times (t^?) 
scale as 

(%V^ /M , f{q)=qD q+1 (41) 

where (.) stands for an ensemble average. However it was found that this relation 
is extremely fragile |S7]. Namely, it holds for channels attached to a typical position 
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Figure 11. Distribution of the delay times V{jw) for the ID Anderson model 
with on-site potential, uniformly distributed between [—0.1; 0.1] and wavenumber 
k = i/n. The dashed line corresponds to 1391 . In the inset we present the same 
data in a log — log plot. The figure is taken from 1841 . 



inside the sample. This excludes the standard scattering set up where the channel is 
attached to the edge of the sample. 

In Fig. 1121 we summarize the results of the investigation undertaken in |87| where 
the analysis was performed for the Power Banded Random Matrix (PBRM) model, 
whose elements are independent random variables i?y with the variance decreasing 
in a power-law fashion: ((Hij) 2 ) = [1 + (\i — j\/b) 2a \~ 1 . For a = 1 this model shows 
critical behavior and the fractal dimensions D q of the eigenfunctions depends on the 
parameter b and can be calculated analytically hi Fig. \12h we report the results 
for the case with a channel attached to the boundary. We see that the numerical data 
deviates from the theoretical predictions for any value of b. Instead, the agreement is 




Figure 12. Scaling properties of inverse moments of delay times for the PBRM 
at criticality 1871 . We report In f(q— 1) as a function of q for a channel attached (a) 
to the boundary and (b) to the bulk of the sample when \ {S) \ ~ (full symbols). 
In (b) we also show ln/(g— 1) for \{S)\ ~ 0.5 (open symbols). The curves are the 
theoretical prediction of Eq. 1411 . The figure is taken from |87) . 

very good for the case where the channel is attached to the bulk of the sample (see 
Fig. 112b ). In the latter case the channel is attached to a representative position in the 
sample. These results will provide a new method for evaluating the fractal dimensions 
D q in microwave and light wave experiments where tw can be extracted even in the 
presence of weak absorption. At the same time the appearance of the anomalous 
scaling of (t^ 9 ) can be used as a criterion for detecting MIT. 
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6. Quasi-periodic systems at criticality 

Periodic and random media cover only the two extremes of the rich spectrum of 
complex systems. Quasi-periodic systems |HH1 IHHI EOl EU EH ESI ESI form 

an intermediate regime and have fascinating properties. In these deterministic non- 
periodic structures translational order is absent. In their one-dimensional tight- 
binding formulation they are described mathematically by the following Hamiltonian 

lpn+1 + tpn-1 + W„1p n = Elp n (42) 

where W n is given by some quasiperiodic sequence. Among the most well studied 
representatives of this class are the Harper model |SS , 93 94 95 and Fibonacci quasi- 
crystals jS^l 191)1 1911 1921 1931 1951 196| . These two systems have been the subject of an 
extensive theoretical and experimental effort in the last twenty years. 

The Harper model is described by the tight-binding Hamiltonian l|42|l with on-site 
potential given by W n = Acos(27r0n). This system effectively describes a particle in 
a two-dimensional periodic potential in a uniform magnetic field with <f> — o?eB/hc 
being the number of flux quanta in a unit cell of area a 2 . When 4> is an irrational 
number, the period of the effective potential W n is incommensurate with the lattice 
period. The states of the corresponding closed system are extended when A < 2, and 
the spectrum consists of bands (ballistic regime). For A > 2 the spectrum is point- 
like and all states are exponentially localized (localized regime) . The most interesting 
case corresponds to A = 2 of the MIT. At this point, the spectrum is a Cantor set with 
fractal (box-counting) dimension Dq < 0.5 [2H1 ESI ES] • The spectral properties of 
the Harper model were recently investigated in microwave experiments |94j . Similar 
theoretical attention was given also to the study of eigenfunctions [SHJ E2j which show 
self-similar fluctuations on all scales. 

The Fibonacci binary quasi-crystal attracted a lot of interest as well. Here the 
potential W n only takes the two values +W and — W arranged in a Fibonacci sequence 
|89) . It was shown [23113 E2] that the spectrum is a Cantor set with zero Lebesgue- 
measure for all W > 0. The first experimental realization of Fibonacci super-lattices 
was reported in ED] while their optical analogues were realized in E3EH1- 

In EH we had presented consequences of the fractal nature of the spectrum in open 
quasi-periodic systems. We considered open systems with one channel (the simplest 
possible scattering problem) and reported the appearance of a new type of resonances 
width and delay time statistics. These distributions show inverse power law behavior 
dictated by the fractal dimension Dq of the spectrum. Specifically, it was found that 
the probability distributions of resonance widths V(T), and of delay times V(t) when 
generated over different energies, behave as 

V(T) =T- a ; a=l + D$ 

V{tw) = ; 7 = 2 - £> E . (43) 

Notice that for Dq — we recover the results (|2(JI39f) of the point-like spectrum 
(Dq = 0) corresponding to a localized system. 

The connection between the exponents a, 7 and the fractal dimension Dq of 
the closed system calls for an argument for its explanation. The following heuristic 
argument [Hi], similar in spirit to the one used in Section 14.11 provides some 
understanding of the power laws (|43|l . We consider successive rational approximants 
4>i = Pi/li of the continued fraction expansion of </>. On a length scale qi the periodicity 
of the potential is not manifest and we may assume that the variance of a wave packet 
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spreads as var(t) ~ t 2D ° [53]. We attach the lead at the end of the segment q^ which 
results in broadening the energy levels by a width T. The maximum time needed 

1/D E 

for a particle to recognize the existence of the leads, is T qi ~ q, . The latter is 
related to the minimum level width T qi ~ 1 /r qi . The number of states living in 
the interval is ~ qi and thus determines the number of states with resonance widths 
r > l/r qi . Thus Pint(r gj ) ~ q% ~ r~ D ° . By repeating the same argument for higher 
approximants 4>i+i — Pi+i/qi+i we conclude that V(T) ~ r~( 1+D ° \ in agreement with 
I l43j l . Furthermore, use of Eq. (|33[1 of Section I5"D results in the distribution H43JI for 
the Wigner delay times. The numerical results (see Figs. 12 and 13) obtained for the 
Harper and the Fibonacci models in [53], verify the validity of the above arguments. 
Nevertheless a rigorous mathematical proof is still lacking. 

7. Conclusions 

In this contribution we summarized the recent activity on the statistical properties of 
resonances and delay times of random/chaotic systems and analyze the deviations from 
the universal RMT predictions due to effects related to Anderson localization, diffusion 
and criticality. We have found that the tails of the resonance width distribution 'P(r) 
reflects the nature of the dynamics associated with the corresponding closed system 
as it is defined by the second moment of a spreading wavepacket (ballistic, diffusive, 
critical, or localized). Instead, the origin, corresponding to small resonance widths, 
is dictated by anomalously localized states. Moreover, we have found that in the 
diffusive regime V(T) is affected by the channel "configuration" (position and relative 
number) as well, in contrast to the localized regime. At MIT the resonance width 
distribution is universal and can be used to formulate a scaling theory of localization. 

Localization phenomena affect also the delay times leaving their traces to the 
the distributions V{t). For scattering systems attached to one open channel, we 
have a very good quantitative understanding of V[t) for all regimes. A general 
expression connects this distribution with the distribution of wavefunction intensities, 
the latter being well studied during the last years. This connection allow us to use the 
experimentally accessible delay times in order to probe properties of wavefunctions, 
like multifractality, which are not easily measured. For many open channels, ample 
of numerical data supported by theoretical arguments allow us to estimate the shape 
of the distribution of delay times and get a qualitative understanding of the traces of 
localization. Nevertheless it remains a challenge to get some quantitative expressions 
as well. 

The last section of this review is devoted to the Id quasi-periodic systems at 
criticality. The corresponding closed systems show fascinating properties like spectral 
and wavefunctions fractality. We reviewed the traces of spectral fractality to the 
distributions of the resonance widths "P(r), and of delay times V(t). Based on 
numerical results and theoretical arguments it is shown that both quantities decay 
algebraically with powers which are related to the fractal (box counting) dimension 
Dq of the spectrum of the closed system. 
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Figure 13. (a)P;„ t (r) of the Harper model (A = 2) for three approximates of 
<T G , <n = °2 = §§§| ; and a 3 = ffff An inverse power law P jnt (r) ~ T 1 "" 

is evident. A least squares fit yields a Ki 1.5 in accordance with Dq ~ 0.5 and 
Eqn. 02). As is seen the lower cutoff of the scaling region decreases for higher 
approximants.(b) Vi n t(r) of the Harper model (A = 2) for three approximants of 
the golden mean a\ = ; <J2 = ®®J 7 ; and 0-3 = j^ 2 ^'^ ■ An inverse power law 
Pintij) ~ t 1-7 is evident. A least squares fit yields 7 S3 1.5 in accordance with 
~ 0.5 and Eq. [Eft. As is seen the upper cutoff of the scaling region increases 
for higher approximants. The figure is taken from 1951 . 
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Figure 14. Power law exponents a, 7 (plotted as a — 1 and 2 — 7) of the 
resonance widths and of the delay time distributions, respectively, as a function 
of the potential strength W for the Fibonacci model. We also plot the fractal 
dimension Dq of the spectrum (the solid line is to guide the eye). Our numerical 
data show that a and 7 are related to the Hausdorff dimension Dq according to 
Eqs. 1431 , The figure is taken from 1951 . 
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